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CHAPTER 1 


INTRODUCTION 


Design of a system to meet a demand or need is a function which must 
be addressed in its entirety. The total system must be examined and the 
interplay of the design variables and system parameters must be intimate- 
ly considered. This systems approach is not the usual case today. The 
system approach is often attempted after the system is ready for delivery 
or for redesign. 

In major systems the piece-meal design procedures used are often 
the cause of inadequate system performance. The systems and units often 
cannot be economically maintained or operated in the real world environ- 
ment. The logistics and support for the system is poor and generally 
determined after the system comes into being. 

Performance, cost, and supportability can be considered key ingre- 
dients in the development of any system. Optimization of these factors 
simultaneously is a difficult task. Trade-offs are always required. 
There exists a need for more interplay between design, system analysis, 
and supportability. The goals established in the areas of performance 
requirements in terms of effectiveness, cost, and supportability must 
initially be addressed in the design process. Design and maintainability 
engineering must translate goals into effective means to make the con- 
cepts viable. This is an interative time consuming process, but impor- 


tant as a first step. 


I 
; 





This study developed a model and algorithm to be utilized as a 
basis for optimal design of a finite population queuing system deployed 
to meet a demand. The sensitivity and constraints of the system para- 


meters and design variables are then investigated with the algorithm. 


Operational System 


The system under study is illustrated in Figure 1, and is designed 
to meet the demand, D. Units in the system are in either of two cate- 
gories; they are available operationally (shaded area) or they are in a 
repair state (in queue or service). In this model the units are all 
repairable when they fail. It is assumed that there are no units beyond 
economical repair or units destroyed. 

All units which fail enter either the intermediate repair level 
(IRL) or the depot repair level (DRL) and will be returned to operational 
use. This model will determine the optimal number of repair channels at 
both the IRL and DRL, the optimum number of units which must be provided 
to meet the demand utilizing the two level maintenance concept, and the 
optimum replacement age. If it is calculated to be more economical to 
have no depot level but only an IRL, the output of the algorithm will 
specify this alternative. The model is flexible in that allowances are 
made to permit a change of the pertinent parameters in order to determine 
the optimum design for specific situations. 

One of the model's primary outputs is the determination of the 
optimum throw-away age. In Figure 1, new units enter the system and 


remain until the nth year. The number of units is assumed to be a 


f 
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constant and the number in each age is also a constant equal to the 
ratio of the total population size and the optimum throw-away age. 

A system approach must be utilized to analyze and integrate all 
the various components and variables involved. In the sensitivity anal- 
ysis of the model, the effects of mean-time-between-failure (MIBF), 
mean-time-to-repair (MTTR), percent allocation for the IRL and DRL, 
depreciation mode, interest rate, and depot MTIR will be investigated. 
In the constraint analysis, shorter desired life for the units and 
population size variations will be analyzed. 

Ihe system approximates the operation of many systems in existance. 
An automobile rental service, numbers of aircraft in a fleet necessary 
to meet operational missions, or industrial equipment required to pro- 
duce specified outputs are a few examples. This model is flexible 
enough to fit many systems with the basic overall characteristics de- 


scribed above. 


Scope and Assumptions 


The model deals with the design variables and system parameters 
which have an impact on the ability of the system to meet the demand. 
However, the system transient phase is not considered, this is a steady 
state model. The portion of the life-cycle to be considered is the 
steady state operation which occurs between the introduction phase and 
the phase-out phase of the system. 

Other assumptions which are made concerning the system are as 


follows: 





1) 


2) 


3) 


4) 


5) 


6) 


The 


not 


The 


but 


All 


number of units deployed to meet the demand is 
constrained. 

population may be made up of units of many ages 
are otherwise homogenous. 


failed units are repairable, no failed units 


will be discarded from the system. 


MTBF and MITR change constantly for a unit, however 


they are assumed to be constant over each yearly 


period, the figures used are yearly averages. 


The 


population is considered finite for calculation 


purposes, and will be held below 50 units due to 


computer limitations. 


Spare units will be in the maintenance system for 


training but will not be considered as a part of the 


total population available to meet the demand. 


system Effectiveness 


The mathematical model formulated in this study has as its basic 


goal the determination of a system configuration which will effectively 


meet the demand. The measure of its ability to meet this demand 


effectively is the cost incurred. The cost equations will be based on 


the annual equivalent cost concepts of engineering economics. 


There are two classes of variables involved in this system, those 


under direct control of the designer and those that are not directly 


under his control. Both types are of primary concern in this development. 





Those which are directly under the designers control are design variables 
whereas those not under his direct control are system parameters. System 
effectiveness depends directly on the values chosen for the design 


variables in the face of system parameters. 


System Parameters 


The eleven variables not directly under the control of the designer 
are (1) initial cost of the units, (2) demand, (3) initial facility 
cost, (4) MTBF, (5) MTTR, (6) percent allocation of failed units to IRL 
and DRL, (7) depreciation method, (8) depot MTTR, (9) interest rate, 
(10) service distribution, and (11) the failure distribution. 

The cost of spare parts, training, documentation, spare units for 
training are also system parameters in this model. The cost of unmet 
demands is a function of the interaction of these parameters and the 
quality and quantity of maintenance and spares. 

The MTBF and MTTR are not directly controlled by the designer, 
however he may decide to have them varied by paying the appropriate 
costs of redesign and engineering to either improve or degrade avail- 
ability. This study does not take into account scheduled maintenance. 
The MTTR will include the transit time and repair time but not the 


waiting time in the queue. 


Design Variables 


The four design variables are (1) throw-away age of the units, (2) 


population size, (3) the number of IRL channels, and (4) the number, of 





DRL channels. The number of units within an age group will be assumed 
equal. For example, a system containing a population of units with a 
maximum age of three years comes into existence by procuring one third 
of the units in each of three consecutive years and replacing one third 
each year thereafter. Under this assumption a population exists which 
will exert a steady load on the maintenance facilities. It will be 
assumed that the system has reached its steady state by some prior pro- 
eurement action. 

The number of service channels to be determined for both IRL and 
DRL is of critical importance. The severity of the failure will determine 
which level of maintenance the unit will enter; the channels at the IRL 
and DRL will have their own repair times. Each channel will have equal 
capability so that the repair time is equivalent regardless of the 
type of failure. The arrivals will be as determined by the proportion 
of failures relegated to each level. The arrivals will enter queues at 
the appropriate facility and be accepted on a first come first served 
basis. The channels required in either IRL or DRL will be parallel 
service channels and they can serve units simultaneously. Demand is not 
a design variable as it is the driving force for the development of the 
system. It is the requirement for which the system is being designed 


to meet. 


Survey of the Literature 


The most significant documents dealing with the procurement and 


acquisition of units or entire systems have been written for the 





Department of Defense. They are "Life Cycle Costing Procurement Guide" 
(4), and "Life Cycle Costing Guide for System Acquisition” (5). The 
former guide is designed to aid in the procurement of the units and 
materials which make up only components of the total system; the latter 
is designed to apply the life cycle cost concept to the entire system. 
These models are applicable to decision making analysis; i.e., comparative 
cost analysis of alternative systems, and analysis of the same system 
under various support or replenishment concepts. These models 

attempt to determine the optimum replacement age for the systems not 
the units in the system. The number of units to be acquired must be 
specified in both of these models. 

There are a number of simulation models available to investigate 
life cycle costs of a system. These models attempt to find the optimum 
number of repair facilities, number of spares required, and personnel 
required to support a specified system. The 'System Cost and Operational 
Resources Evaluation (SCORE)" model developed by the Naval Air Develop- 
ment Center (12), is an example of such a model. It estimates the life 
cycle costs of research and development, test and evaluation, investment 
and operations for fifteen years for the various component elements, 
and aggregates these into a total cost estimate for the system. The 
models mentioned above require large amounts of data to adequately per- 
form their analysis and require extremely long computer runs. 

Much of the effort in the area of system costs has been directed 
to defining the reliability and maintainability aspects of the systems. 


Armstadter (1), and Gnedenko (8) do mathematical analysis of systems 





to attain the probable best combinations of reliability and maintain- 
ability but do not consider any cost functions and their interrelation- 
ships to these parameters and effects on the system performance. They 
define MIBF as a measure of reliability and MITR as a measure of main- 
tainability. Seiler (15) does an analysis of a system assuming that the 
probability functions of cost are triangular in shape and that there is 
only one optimal cost with all other costs linearly distributed on 
either side of the optimum. Briggs (3) did work in determining the 
reliability of systems and its effect of cost. His effort is in areas 
which are generally post design phases, after the system has been in 
operation for a period of time. 

MILHANDBOOK 472 (11) defines the various distributions of failure 
and repair times and rates. Specific statistical distributions are 
listed which correspond to those used in the model. The basis for the 
queuing equations used in this model is similar to that described and 
utilized by Jelen (9) and also those problems in the area of queuing 
as discussed by Kendall (10). 

Wortham (17) bases his analysis of system cost on the initial cost 
of the system and support costs. In his analysis, the availability of 
the individual units was considered as not critical when the value was 
between 80 and 85 percent because then there would exist ample time to 
perform maintenance and support functions. However if the availability 
was necessarily between 95 and 97 percent, then the design reliability 


and maintainability were considered very critical. 





1.0 


In Goldman (7), the use of a particular definition as to what are 
support costs can change the various component life cycle cost calcula- 
tions considerably. Many large systems support costs often exceed the 
cost of the equipment by more than ten times, even 100 times is not rare. 
Also, Brandt (2) developed a model which is utilized to predict the 
down-time of a unit as a function of its utilization, and whether it is 
utilized fully or partially when operated. 

In summary, a review of the available literature on subjects 
related to this thesis does not reveal any method of optimizing all 
the design variables for a general operating system as addressed herein. 
Further, none of the analyses consider the level of maintenance concept 


or the impact of the various variables on system effectiveness. 





CHAPTER 2 


MATHEMATICAL MODEL AND COMPUTER ALGORITHM 


In attaining the primary goal of this study it was necessary to 


develop a model and an associated algorithm which would permit determi- 


nation of the throw-away age, population size, and number of service 


channels in IRL and DRL to meet demand in an optimum manner. The model 


is cost based; it incorporates the annual equivalent cost (AEC) of the 


maintenance facilities and the annual equivalent cost of the units plus 


the annual cost for non-fulfillment of demands. 


The following assumptions are made in the development of the 


mathematical model: 


1) 


2) 


25 


4) 


5) 


Inter-arrival times are exponentially distributed at 
both repair levels. 

The repair times are exponentially distributed at 
both repair levels. 

The inter-arrival times are statistically independent 
of the arrival times, at both repair levels. 

Ihe service channels, at both IRL and DRL are parallel 
and each is capable of performing the maintenance 
activities associated with its respective area of 
responsibility. 

The population size will always be larger than the 


number of service channels. 


111. 
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6) Each channel performs service on one unit at a time. 

7) Units are serviced in the repair channels in a 
first come first served mode. 

8) Units completing repair return to operation with the 
same operational, reliability and maintainability 
characteristics established for its associated age. 

9) The cost of the repair facility varies linearly with 
the number of channels. The initial investment cost 
in land and associated improvements is constant for 
up to three channels; the base cost of the facilities 
is increased for each additional increment of three 
channels, this pertains to both levels of maintenance. 

10) The out-of-operation cost is a linear function of the 
number of units not available to meet demand. 

11) The administrative, labor, and operating expenses for 
the IRL and DRL facilities increase linearly with the 
number of channels being required. 

12) The cost of training, initial spares, and recurring 
support are included as proportions of the total annual 


equivalent cost (TUAEC) and the IRL facility costs. 


Population Annual Equivalent Cost 


To determine the cost of ownership of the population of units on 
an annual equivalent cost basis, the fundamental equations of engi- 


neering economy are used with the notation modified for this thesis (16). 





5 


TUAEC = Population annual equivalent cost 


D = First or acquisition cost of the unit 

SV = Estimated salvage value of the unit 

BT = Book value at end of year AGE 

AGE = Year of oldest throw-away age, AGE < NP 

NP = Estimated life of the unit 

N = Number of units in the population 

፲ = Interest rate 

TUAEC = {(P-BT) ian _ + BT(i) $ N 
(1+1) -1 


The book value of the units at the time of their disposal will 
depend upon the above equation and the method of depreciation used. 
The methods of calculating the depreciation of an item over time are 
given in the following sections. 

Straight line method. The straight line method of depreciation 


assumes that the value of an asset decreases at a constant rate. 


BT = P - AGE ES 


Sum of digits method. The sum of the years digits depreciation 
method assumes that the value of an asset decreases at a decreasing 


rate. 


NP-AGE NP-AGE+1 
“5.3 Ve 


BLEU DEL EE * SV r 
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Sinking fund method. The sinking fund method of depreciation 


assumes that the value of an asset decreases at an increasing rate. 


i Mo MM 
0+0 1 98٦ 


BT = P - (P-SV) 
Zero book value. Although not specifically mentioned in the 
references, in some circumstances it is desirable to set a zero value 
on an asset. This will take into account the situation where units 
are retired (at any age) with no appreciable value. From the above 
equations, the total annual equivalent cost of all units can be 


calculated knowing the age of the oldest units. 


Repair Facility Annual Equivalent Cost 


The annual equivalent cost for both IRL and DRL facilities must 
be considered. The basic formula and the notation to be used are as 
follows: 


eles) Ee 


FAEC = (PFAC-FFAC) 
NEE 


+- FFAC (i) 
where, 
FAEC = facility annual equivalent cost per channel 
PFAC = first cost of the facility per channel 
FFAC = estimated salvage value of the facility per channel 
LFS = estimated life of the facility 


i = interest rate r 





15 


This equation is similar to the one derived previously for TUAEC. This 
cost is a linear function of the number of repair channels and has a 
non-zero, positive value when the number of repair channels are cal- 
culated. The value of the depot facility annual equivalent cost (FAECD) 
for the situation where there are zero depot channels is zero, there 
will be no investment of any type at that level. However, the cost of 
the IRL will be increased in the algorithm to take into account the 
increased repair requirements at that level as its maintenance respon- 
sibility would be much greater. 

In the calculation of the facility cost the various types of 
equipment and buildings must be taken into consideration, these 
different values and different life expectancies are handled as a 
summation for each level of repair. Each facility can handle only a 
given number of channels. A factor is added to the algorithm which 
will increase cost as the number of channels is increased. This cost 


is a function of the number of channels for both types of facilities. 


Let 
FAECN = total facility annual equivalent cost 
DY = yearly administrative, manpower, and overhead 
costs per channel 
CC = number of IRL channels; CC > 0 
(for DRL, DC = number of DRL channels; 0 < DC < 6) 
ALEO = first cost of establishing facilities 
Then, e 


FAECN = CC(FAEC + DY) + ALEO 
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This equation is utilized with two sets of notation in the model 
to facilitate handling of TRL and DRL. The notation for FAEC which 


contains a D, pertains to the depot facilities. 


Out-of-Operation Cost 


When units are not available to meet a specified demand some cost 
is incurred. The number of units available at any given time is a 
random variable, but is of unknown distribution. The out-of-operation 
cost is therefore a function of the demand and the number of units 
available. This number is a function of the failure rate, the time to 
repair, the population size, the number of unit age groups, the expect- 
ed ratio of major and minor failures, and the number of IRL and DRL 
facilities being provided. 

The first parameters which need to be determined are the MTBF and 
MTTR values at both levels of maintenance, for the populations at each 
level. It is assumed that the value of these parameters is known or 
can be estimated initially. They are input system parameters and 
extend over the life of the units. These inputs are on a yearly basis. 

Since it has been assumed that the number of failures per unit 
time and the repairs per unit time are Poisson distributed, a theorem 
by Parzen (13) can be used for the cumulative sums for the various 
unit ages. The theorem states that if S, = Xy FXy ho KL is the sum of 
n independent random variables, if k = 1,...,n; Xx is Poisson distri- 
buted with parameter A rot Ç Also the mean of S, is equal to the 


1 


sum of the individual means. e 
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Let 
SEH = mean time between failures for the jth unit age 
nun, = mean time to repair value for the jth unit age 
n = number in the unit age 
N = population size 


Assuming each throw-away age contains the same number of units, 


the calculations necessary to find the MIBF and MTTR are as follows: 


n wm, (S 
MTBF, = ) dd (2:53) 
JE N 
j=1 
n NIR. S 
MITR, = (232) 
j=l 


Once these values are determined, the expected number of units 
operational can be determined by using the Finite Queuing Tables of 
Peck and Hazlewood (14) or computed using the equations developed in 
Appendix A. The results derived in Appendix A are used for this 


algorithm and the notation given below appears in the model. 


H 
wae MTBF 
اط‎ g MTTR 
H 
p. E 
y 
W = expected time a unit waits in a queue ' 
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n 
= ce( f 
SE =ርርተ1 N! SZ ን 
የበኩ ፡።ይ o e 
° cc (N-n)! cc! 
W = BECC D n 
q cc-l N! : n N! ce (ራ) 
En) en) mr -- 
n=0 nI (N-n)! SE (N-n)! cc! 
(2.3) 
L = expected number of units in a repair level (queue 
plus repair channel) 
: _ 4፲/ህተጊ/ለን ፒጅ zellen 
Ho 1s (1/ህ+1/እ4፪.ን S O ENS 7 


The expected number of units available is determined using the 


following equation: 


TJ = XN - (L + LL) 


where, 
TJ = expected number of units available 
XN = total unit population 
L = expected number of units in IRL 
LL = expected number of units in DRL 


This number is utilized in calculating the cost incurred for not meeting 
a demand. Since the number of units available is a random variable, it 
has some type of distribution (which is not determined herein). It 

is assumed that there is a degree of variability in the number of units 
available which could alter the cost of not meeting the demand. For 


simplicity this cost will be assumed to be a linear function of the 
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number of units demanded and not available. 


Let 
C. = cost of unmet unit demands per day 
Z, = cost of unmet demand per unit time period 
D = expected demand 
Q = degree of variability o£ the number of units 


available at any given time 


(1J-D), Q> (TJ - 5)‏ ےہ کٹ 


0 O CGD) 


The system cost includes the training of personnel to operate and 
maintain the units which are procured for the system. Seven percent 
of the total unit and facility annual equivalent cost is added to the 
system cost for training and initial support. The spare parts and 
piece parts are a recurring expense, twenty percent of the TAEC of the 
units and facilities is used in this model as their cost estimate. 
The grand total annual equivalent cost is the sum of all costs pre- 


viously specified and is expressed as follows: 
TIAECN = TIAEST + FFAEST + STRN + RSUPN + WCN 


where 


TTAECN 


grand total annual equivalent cost 


TTAEST 


total annual equivalent cost of intermediate 


plus depot level units 
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FFAEST = total annual equivalent cost of the intermediate 
plus depot facilities 
STRN = training and initial support cost 


RSUPN = recurring support cost 


WCN = cost due to not meeting a demand 


The Computer Algorithm 


The program coding which was developed for this algorithm is 
exhibited in Appendix B.1. The program is made up of a main routine 
and three subroutines. Each of these areas will be described in the 
following sections. 

Main program. The primary objective of this program is to minimize 
(optimize) the TAEC. The goal of the main program is to execute a 
search technique which will optimize the total cost function over the 
range of variables delineated by the designer, and to act as the handler 
of the other routines. The total system cost function is unimodal. 
There is only one true minimum for any given set of system parameters. 
The number of variables over which the search must progress is four: 
the number of IRL channels, the number of DRL channels, total popu- 
lation size, and the throw-away age of the units. 

To commence the search it is necessary to determine the service 
factors to be used in the IRL and DRL. This equation was originally 


developed and used by Peck and Hazlewood (14). 
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Let 


MuR E 
X = Service factor (IRL) XMTTR + XMTBF 
XD = Service factor (DRL) PD 


DMTTR + DMTBF 


Using the number of channels in IRL as the starting function for 
the search routine, the number of channels from which to commence the 


search is calculated as follows: 
CC = (N - X) + á [1 * ፤(1-3)1”"* 


The search is run through the range of IRL channels looking for 
the minimum TAEC. During this search the number of DRL channels is 
maintaineá at a constant value. After each minimum TAEC is determined 
for a population size, the DRL channels available is incremented. A 
run is made with each succeeding DRL value. Only the minimum set of 
TAEC data for each combination of IRL and DRL channels is stored in 
the computer. This minimum can be visualized as a point on a plane in 
three dimensional space, the plane being for a given population size 
and throw-away age. For each throw-away age considered there will be 
a set of outputs; these must be visually scanned for the minimum TAEC 
for that throw-away age. When all population sizes and throw-away ages 
which are being considered have been searched, the minima must be 
checked again to identify the fourth minimum; this is the ultimate 


objective of the main program, the global minimum TAEC. The results 


7 
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indicate the number of IRL AND DRL channels, and the associated cost 
of that population size and throw-away age. 

Subroutine RAIDS. This subroutine determines the rate of arrivals 
and service at IRL and DRL. It provides the aggregate XMTBF and 
XMTTR, DMTBF and DMTTR values for IRL and DRL for each population size 
and throw-away age combination each time they are required in the com- 
putations. The equations in (2.1) and (2.2) are utilized in the 
execution of this subroutine. This routine includes the calculation 
of the inherent availability of the units in the intermediate level. 
This is the designed availability (not the actual availability) which 
is determined after the system has been in use for some period of time. 
The availability is contained in the printout for each population size. 


It is computed as follows: 


E XMTBF 
XAL = TBF + XMTIR 
XAT = Inherent availability 


Subroutine UNAVAIL. This subroutine determines the number of 
units available from a given population size utilizing the appropriate 
values determined in RAIDS. The important result of the calculations 
using equations (2.3) and (2.4) is the expected number of units not 
available. This subroutine is limited in the population size it can 
accommodate due to the size of the factorials which must be calculated; 
55 factorial is the limit for the computer used with this model. A 


limit may also be reached if the desired availability of the units is 





A 
ከ. 
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high. Generally the program will function adequately when the ratio 
of MITR to MITR plus MTBF is greater than two tenths. The maximum 
population size which can be run on the IBM-370 Computer System irre- 


spective of the ratio described above is 43 units. 


Subroutine INVEST, All cost factors mentioned previously are 
included in the overall system cost and are calculated and summed in 
this routine. It includes the MOD designation for the book value 
calculations. This routine calculates all the costs for the IRL and 
DRL facilities which includes an escalation factor for the increased 
number of channels in either facility. The TUAEC for the units in both 
levels of maintenance, training and initial support, and recurring 
support costs are also computed in this routine. All costs are totaled 
and the result is the TAEC for the particular combination of para- 
meters being investigated. The basic time unit used throughout all 


the routines is one day in a 365 day year. 


EN 


l ۲ 





CHAPTER 3 


STANDARD SYSTEM DESIGN EXAMPLE 


This chapter presents a standard system design formulated to 
utilize the mathematical model and exercise the optimization algorithm 
developed in the previous chapter. The standard system parameter in- 
puts must be obtained by the designer and represent a system with which 
he is concerned. He will receive from the programmed algorithm the 
optimal combination of IRL and DRL channels, the minimum TAEC, and when 
to throw-away the units for the most effective and efficient system 


٥٦٦۳67٦12780610۰ 


System Parameter Input Values 


Appendix B.2 gives the standard system design input values. The 
demand is 20 units per day. The maximum age of the units in the system 
is 15 years. The depreciation code is two (straight line method). The 
number of units to be repaired at each level is allocated 75 percent to 
IRL and 25 percent to DRL. 

The depot MITR is 2.5 days and is constant over the life of the 
units. All units, when sent to depot, will be overhauled and returned 
to service with the same performance and maintainability character- 
istics as their original age group. The depot MTBF will be the same 
as the IRL MTBF; all units will have the same reliability. 

Ihe severity of the failures is not a function of the MTBF but 


is a random factor in the unit design which is unpredictable. The 
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value of MTBF and MITR for the IRL are given next. A value for each 
year is estimated and is the average value for the units in that throw- 
away age, i.e., in the second year the MIBF is 6.5 days and the MTTR 

is 1.3 days. 

The MTBF of the units in this system are assumed to have the "bath 
tub" failure hazard characteristic. During the first few years the 
MTBF will improve and reach its designed value. It will slowly degrade 
as the units age. As the units near the end of their projected life 
the failure rate will increase rapidly. 

The MITR will also vary. At initial introduction the time to 
repair will be longer than desired for the units and will reach its 
design value about the same time that the MTBF achieves its design 
value. The MTTR will increase as the units age due to in-depth repairs 
becoming more prevalent. 

The cost data used in this example will be as follows: The 
initial unit cost is $25,000, with a maximum life of 15 years. The 
facility costs are separated into three categories for both the IRL 
and DRL. The first of each set is for buildings and related equipment, 
with the IRL expected to last 40 years and the DRL 30 years. The 
large repair equipment, test benches, and test equipment is the second 
set and will be usable for 10 and 8 years respectively. The third set 
is for tools, test kits, and small support equipment and will last 5 
years at both levels. All facility costs are listed in Appendix B.2 
and are calculated in the subroutine INVEST on a per channel basis. 


The interest rate on all investments is 10 percent. e 
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The expected yearly administrative labor and overhead costs are 
$25,000 and $62,500 respectively at IRL and DRL. The cost of land, 
utilities, etc. are included as inputs and are computed in the algorithm 
as a step function dependent on multiples of three repair channels. 

The initially assumed cost is $25,000 and $75,000 respectively. 

The cost for not meeting a demand is $200 per unit per day. The 
units available will vary by a factor of four when there is no DRL 
and by a factor of 2.5 when both levels are available. If the number 
of available units is equal to or greater than 22.5 there will be no 
cost incurred, however should the number of units available be less 


than 22.5, a cost will be incurred. 


Design Variable Output Results 


Appendix B.3 gives the results using the algorithm and standard 
system design example input. The printout shows all the costs which 
make up the minimum cost combination of IRL and DRL channels for a 
given throw-away age and population size. To allow for better under- 
Standing of the output of the model, the seventh throw-away age 
listing and a population of 34 units will be discussed in detail, then 
the entire throw-away age range is analyzed. 

The algorithm starts at throw-away age one and proceeds to the 
maximum age designated in the input of 15 years. The unit inherent 
availability for this case is 86.21 percent. The number of IRL 
channels is 5 and the DRL channels is 1. The queuing routine described 


in Chapter 2 is used to determine the number of units that will not 
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be available. The channels are searched until a minimum TAEC is found. 
For this problem the cost computed for IRL is as follows:  TUAEC is 
$106,083.81, and the FAEC is $176,957.88. For the DRL, the TUAEC is 
$35,361.28 and the FAEC is $15,006.44. The training and initial support 
cost is $19,812.91, while the recurring support cost is $56,608.33. The 
cost for not meeting the demand is $0.00. Also listed are the number 

of units allocated to IRL (26) and DRL (9). The minimum total annual 
equivalent cost is $544,830.56. 

A search of the various combinations of IRL and DRL in a throw- 
away age is executed and the optimum for a population size is printed 
out. In each population the number of IRL channels searched will 
normally be between 1 and 12, while the number of DRL channels can 
vary from 0 through 6. There must be at least one IRL channel; main- 
tenance must be performed at that level but there need not be a depot 
level of maintenance. 

Once a printout is available the minimum cost in each throw-away 
age can be determined. The minimum of these minima is the global 
minimum. In the printout of the standard example, the throw-away age 
containing the minimum is 10. The total population size is 39, and 
the TAEC is $519,324.81. 

Figure 2 shows the minimum TAEC versus throw-away age, this plot 
shows which throw-away age has the minimum cost. The associated print- 
out must be used to determine the system configuration which produces 


this minimum cost. 
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Figures 3 and 4 are graphic presentations of total population size 
versus the minimum total annual equivalent cost. Each throw-away age 
is represented and the global minimum can be seen in Figure 4. These 
figures are developed from the standard run whose output is included in 
Appendix B.3. Two figures are used to facilitate reading with two 


different vertical scales to display the results. 


Discussion 


Analysis of the output in Appendix B.3 indicates that there is 
considerable information available to the designer beyond that already 
discussed. Should he have the option to use a range above the minimum 
TAEC for the system to be procured, there are several system con- 
figurations which now become acceptable. In this problem if the limit 
was plus 2 percent of the minimum cost, he would then look for all com- 
binations with a cost less than $530,400. Scanning the output data it 
is seen that there are 14 options within the 2 percent range limit. 

Another area the designer may wish to investigate is when there 
is a maximum amount of financial resources available or a limit on the 
acceptable annual equivalent cost. He can quickly ascertain from the 
Output if there is an acceptable system configuration. He may need to 
rerun the program varying the system parameters until he attains an 
acceptable solution or decides that he cannot meet the demand with an 


acceptable system configuration. 
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CHAPTER 4 


ALGORITHM SENSITIVITY ANALYSIS 


This chapter presents an analysis of the mathematical model and 
the optimization algorithm developed in Chapter 2. Several system 
parameters will be varied and their effects on cost will be investigated. 
This analysis will allow the identification of critical parameters. 

The parameters to be tested will have their values specified by 
the designer. These parameters will be assumed to be logical estimates 
and will be substituted into the standard input set. The original 
system design example output will be used as the basis for all compari- 
sons. The variable values used in making the test runs will be listed 
with each test discussion. The lines in the standard input set, 
Appendix B.2, will be numbered from top to bottom, one through 15, for 


purposes of reference in this chapter. 


Reliability and Maintainability Variation 


The item to be investigated first is the determination of the 
system sensitivity to changes in MIBF and MTTR. These changed factors 
result in either the increasing or decreasing of unit availability. 
Availability is used to realign the unit cost. The resultant changes 
to be incorporated in the new input set are given below. 

It is assumed that for each one percent increase in availability 
the initial unit cost will increase by $3,000. The cost will be de- 


creased for a drop in availability at the rate of $500 per one percent 
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change. The following is a listing of the system input data used in the 


test runs. 


First Test Run 


Line 4 MTTR: same as standard run 

EX BE (days): 9.1: 10.5; 12.05; 12.0; 12.0; 11.5; 
۳0۰۰۱٣٦۰ ۱100۰95 9.0 8.5: 0,0; 7.0;. 6.0; 

Line 6 Unit cost input: $40,000; $4,000 salvage value; 15 year 


life; Interest rate 10%. 


Second Test Run 


Line 4 MTTR: same as standard run 
mne 5 MIBE (days): 4.0: 5.0; 6.05 6.0: 6.0; 6.5; 6.25; 
OPO O O ba oe 4 Os), EE 
Line 6 Unit cost input: $23,100; $2,310 salvage value; 15 year 


life; Interest rate 10%. 


Third Test Run 


e ETIR da yS datos 155: 1255 T 25 01.25; 1.255 153253 
726 ep u455ure5:001:6; T7: 1:8: 1.9: 2.05 

Line 5 MIBF: same as second test run 

Line 6 Unit cost input: $21,800; $2,180 salvage value; 15 year 


life; Interest rate 10%. 
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Fourth Test Run 


EUST MITR (days): 1.3; 1.2; 1.1; 1.0; 0.9; 1.0; 1.1; 
icm o1 3sclo3oe S4: 1:5: 1.05. 1.9. 
ieee > MIBF (days): 5.5; 6.9; 8.6; 8.6; 8.6; 8,0; 7.5; 
I, 10: 0.5: 0.5: 6.0: 6.0; 5.57 4.5; 
Line 6 Unit cost input: $29,500; $2,950 salvage value; 15 year 


life; Interest rate 10%. 


As indicated in Figure 5, it is more economical to procure and 
deploy a system with the standard system design parameters than any of 
those tested. It is more cost effective to increase the unit avail- 
ability, i.e., the ratio of MTBF and MTBF plus MITR, than to decrease 
it. A negative increment of the availability resulted in a higher 
TAEC than the same increment in the positive direction from the 


standard run availability. In Figure 6 this fact is also shown. 


Percent Population Assigned to IRL and DRL 


In this test the percent of the total population of units in tne 
system allocated to IRL and DRL is varied. The interval is 65 percent 
to 85 percent inclusive, in five percent steps. This allows a repre- 


sentative cross-section of allocation estimates the designer might use. 


First Test Run 


Line 2 IRL percent = 65 
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Second Test Run 


Line 2 IRL percent = 70 


Third Test Run 


Line 2 IRL percent = 80 
Fourth Test Run 
Line 2 IRL percent = 85 


As depicted in Figure 7, the plot of the TAEC versus throw-away 
age indicates that there is very little variation in cost when compared 
to the standard percentage, 75. This analysis shows that the population 
allocation to IRL and DRL could vary between 65 percent and 85 percent 


with only minor changes in the minimum TAEC. 


Depreciation Method 


Ihe next parameter modification is the method of depreciation 
utilized in determining the book value of the units. The standard 
design example run was made using the straight line method. In this 
analysis the other three methods available in the model will be examined. 
First, the sum of years digits; second, the sinking fund method; and 


thirdly the book value equal to zero are investigated. 


First Test Run 


Line 1 Depreciation code = 3 ' 
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Second Test Run 


Line 1 Depreciation code = 4 
Third Test Run 
Line 1 Depreciation code = 5 


In Figure 8 it is readily apparent that the two modes, book value 
equal zero and the sum of years digits method are less costly during 
the initial ownership years. However, they both approach the straight 
line method when the throw-away age approaches the maximum life. The 
sinking fund method is very much more expensive in the initial owner- 
ship years. It too approaches the other methods as the units become 
older. They all are equal at the end of the life cycle of 15 years 


as they are expected to be in this model. 


Interest Rate Variation 


The TAEC will be varied because of incremental changes in the 
interest rate. The effects of these changes will be investigated in 


this section. 


First Test Run 


Lines 6 through 12 Interest rate = 8% 
Second Test Run 
Lines 6 through 12 Interest rate = 9% í 
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Third Test Run 
Lines 6 through 12 Interest rate = 11% 


Figure 9 shows the result of the new input parameters and the 
computations for the TAEC. The variation is equally spaced from the 
standard run cost for a particular throw-away age. This should be 
expected. The difference at the ten year minimum cost level for the 
standard run and the two adjacent percentages is approximately $12,000 


above and below the standard run TAEC. 


Depot Level Repair Time Variation 


In this, the final sensitivity analysis test to be run, the depot 
repair time is varied. This test will identify if there is a limit 
to the time out-of-service due to a major failure which results in the 
unit being sent to the depot facility for repair, and its effect on the 


TAEC. 

First Test Run 

Line 3 DMTTR (days): 5.5 
Second Test Run 

Line 3 DMTTR (days): - 12.5 


The results presented in Figure 10 shows that there is an apprecia- 
ble increase in system TAEC when the mean-time-to repair is increased 


by 800 percent. However, slight changes in DMTTR, up to about 5.5 days, 
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cause only moderate increases in TAEC, this is not considered an ex- 
cessive variation. This fact is of importance to the designer as it 
allows for some latitude in the depot facility design and location as 
well as the amount of test equipment and level of spares that must be 


on-hand for effective and efficient unit repairs. 


Summar y 


A summary of the results of the sensitivity tests is given in 
Table I. The table lists the percent change in the parameter varied, 
with respect to the standard run parameter, and the resultant change 


in the TAEC, compared with the minimum TAEC. 
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TABLE I. Percent Variability of TAEC in Sensitivity Analysis 


Throw-away Age = 10 % (STANDARD) 












Reliability and Maintainability Std. = 85.9% 





Run No. 1 
Run No. 2 
Run No. 3 
Run No. 4 












+07 


Percent Population IRL/DRL Std. 75% 















Run No. 1 
Run No. 2 
Run No. 3 
Run No. 4 





Std. = 2 


Std. = 10% 


Method of Depreciation 














Run No. 1 
Fun No. 2 
Run No, 3 





Interest Rate Variable 





Run No. 1 
Rün NO., 2 
Run No. 3 


DMTTR Variable Std. + 2.5 


Rum NO. 1 +220 
Run No, 2 +800 





4 (TAEC) 


$519324.81 


$519324.81 


+209 
12.09 
05 
-0.79 


5519324 .81 


-4.76 
TOOL 
-2.66 


$519324.81 


-4.14 
-2.09 
72.10 


5519324 .81 


+2.43 
+5.12 





CHAPTER 5 


ALGORITHM CONSTRAINT ANALYSIS 


This chapter will extend the analysis of the model to investigate 
the application of constraints on the design variables to determine 
their effects on system effectiveness. The two variables that will be 
considered are the throw-away age and the total population size. The 
format for the test run discussions will be the same as that used in 


Chapter 4. 


Throw-away Age Constrained 


The throw-away age can be restricted by the designer or by manage- 
ment to only have units in the system newer than a prescribed age. In 
this case he would want to evaluate the effect of such a decision on 
the TAEC. This investigation requires four of the data input lines 


be modified. They are as listed below. 


First Test Run 


Line 1 NNAGE = A 
Eine 4 MIIR (days): 50252 120: 1.25; 1.9. 
Eine 5 MIBE (days). 5.0: 8.0; 6.5; 4.0. 


Line 6 modify life of units to 4 years. 


Second Test Run 


Line 1 NNAGE = 8 I 
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né 4 MIIR (days): 1.5; 1.1; 0.9; 1.2; 1.3; 1.4; 1.6: 1.9. 
ቪር SU MIBE (days): 5.0; 8.0; 8.0; 7.0; 6.5; 6.0; 5.5; 4.0. 


Line 6 modify life of units to 8 years 


Third Test Run 


Line 1 NNAGE = 11 

ee 4 MITR (days): 1.5; 1.3; 1.0; 0.9; 1.2; 
i235 415401555 los 1.75: 1.9. 

ne o MIBF (days): 5.0; 6.5. 8.0; 8.0; 7.0; 
03952.0202.5293.9:953 9.0; 4.0, 


Line 6 modify life of units to ll years 


In Figure 11 the program output is plotted as previously; TAEC 
versus Throw-away Age. The results are very definite, the cost incurred 
by limiting the unit life is very high. The TAEC does follow the ex- 
pected trend as it will come very near the standard design run as it 
approaches 15 years. At eleven years, with a throw-away age of 10, the 


difference in TAEC is plus $50,000. 


Total Population Constrained 


There may also be some physical constraint placed on the scope of 
the system. The space available in which to store the units could be 
a very serious problem. In this case the designer would have to limit 
the number of units to be procured and ascertain what effect this has 


on the TAEC. In this problem two lines of input data need to be modified. 


7 
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First Test Run 


fane | XN = 35 


Line 2 ENDXN = 26 


Second Test Run 


Line 1 XN = 30 


Line 2 ENDXN = 21 


Figure 12 shows the expected effect of a constraint and also an un- 
expected effect. In the case of the second test run, an increase in the 
system cost is exhibited due to the population size constraint. However, 
the first test run shows the effect of rounding the factorial and com- 
binatorial calculations. An exhaustive analysis of the terms in the 
equations described in Chapter 2, and contained in subroutine UNAVAIL, 
lead to consideration of the terms invoiving factorials and combinator- 
ials. These terms are rounded amd as a result the accuracy of the 
equations is affected. These calculations in turn result in the varia” 
tion of the number of units available. This number is relatively small 
but has a considerable effect on the cost of not meeting demand. As 
an example, a change of one-tenth in the number of units available per 
day below the established minimum required in the algorithm will re- 
sult in a cost of $7,300 per year. In the case analyzed, the unit 
availability changed sufficiently to incur a cost on the system in 
the unconstrained case but not in the constrained thereby giving the 


false indication of a lower system cost. 2 
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summary 


To summarize, the constraint analysis shows that the effect of con- 
straints on system cost is more pronounced and they have greater influence 


on the TAEC than system parameter variations. 
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CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 
As a result of this study it is concluded that a multivariant 
system for a finite population with multiple queue channels, a two 
level maintenance capability, and optimization of the system to meet 
a demand at minimum total system cost can be achieved. In Chapter 3 
it was shown that utilization of the designers parameters does lead 
to the determination of a global minimum. The optimum throw-away age, 
the number of intermediate and depot level facilities, the optimum total 
population size, and the minimum total annual equivalent cost for that 
particular system configuration can be specifically determined coinci- 
dentally. 
From the results of the sensitivity and constraint analyses in 
Chapters 4 and 5 respectively, it can be concluded that: 
1) The MTBF and MTTR values utilized by the system designer 
has considerable effect on the overall TAEC. 
2) The optimal inherent availability of the units to be 
procured should be approximately 86 percent to achieve 
the most effective system for the case studied. 
3) The percent allocation of the failed units to the 
intermediate or depot level of repair has little 


effect on the total system cost in the range analyzed. 
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The depot MTTR has little effect on the TAEC when the 
increase in MITR is less than 200 percent. Large changes 
in MTTR do increase the cost by more than five percent. 
The interest rate and the depreciation mode used in the 
model have predictable results. All the results were 
consistant with the formulas incorporated in the model. 
Limiting the throw-away age of the units will raise the 
TAEC of the system if the limit is low in comparison to 
the unit's original design life. 

The constraining of the total population size, combined 
with the system parameters estimated and utilized in this 
model was such that the effects of computer word size and 
rounding were dramatically exhibited. This is the result 
of fortunate selection of system parameters and varlables. 
The second test run exhibits the normally anticipated 


results. 


In this study, the concept of designing the total system from the 


beginning is considered to be the critical factor in achieving the 


optimal system design. The demand is the starting point in arriving 


at the specification of a system which will optimally meet the demand. 


The total population, the types and number of repair facilities, the 


reliability, maintainability, costs of unmet demands, initial cost of 


units and facilities, and the necessary support and logistics were 


all addressed in this study. 
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The mode of transformation of all the requirements is the model 
and its associated algorithm. This is but a small part of the planning 
and estimating required to bring a system into being efficiently and 
effectively. The designer can now write the specification and allocate 
the requirements to the various engineering and manufacturing specialists 
for their use in the detailed engineering and production required to 


bring the system into being. 


Recommendations 

There are several areas addressed in this study which should be 
analyzed further. They generally pertain to the basic system design. 
A few of these are as follows: 

1) Further refinement of the model is needed in the areas of 
support and personnel costs. The estimates made are 
percentages representative of systems in use today. More 
exact methods should be sought to delineate these costs 
more accurately. 

2) An analysis of the model and necessary modifications to 
the algorithm should be made to enable more detailed 
investigation of variations in demand and constraints on 
the number of IRL and DRL channels and their effect 
on the system total cost. 

3) A study should be made to identify the parameters and 
variables effected when the total system cost is limited 
to some maximum amount. This will require a substantial 


Ld 


change in the search routine. 
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An analysis of the out-of-operation units should be 
made to establish the distribution of the units down. 
This will facilitate establishing a better equation 
for the unmet demand cost. 

Continued efforts should be made to determine the 
critical system parameters early in the system design. 
The algorithm and model can be further refined, by 


updating equations and reducing the number of assumptions. 
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Development of the Finite Population 


Queuing Equations 
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The following mathematical development is contained in Chapter 111 
of a manuscript entitled, "Analysis of Queuing Systems" by G. K. Bennett, 
J. W. Schmidt, and J. A. White of the Department of Industrial Engineering 
and Operations Research at Virginia Polytechnic Institute and State 
University. 

This development utilizes the classification technique introduced 
by Kendal1(10). The basic queue developed for use in this thesis is a 
Finite Queue, (M/M/c):(GD/N/N) is the type which is developed in this 
Appendix. In order of occurrence, the letters used above are interpreted 
as follows: 

M - denotes Poisson arrivals to the repair facility 

M - denotes Poisson repair times 

c - denotes the number of repair channels 

GD - indicates the units are served according to any general 
queuing discipline 
N - is the maximum number of units allowed in the repair 
system (in repair and waiting) 
and N - is the finite population size. 
The probability mass function for the number of units in the system 


is as follows: 


T 7 ነ . 
À = MTBF " expected number of failures per unit time 
። m = expected numbe f repai er unit time 
u MTTR xp umber of repairs per 
À . . # 
ዐ= ae = traffic intensity, system utilization factor, or 


system load factor. 
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n n 
E Ea n=0,1,...,c+l 
c n 
Ny nic 
E ۶ھ“‎ S n=c,ctl,...,N 
0 ; n> N 


Solving for Po we have: 


c-1 N ۰ N N MS n “3 
o š 3 e (ep) di de (2) e 


To determine the expected number of units in repair, L, and the expected 
time a unit waits in the queue, E let the following equation be used; 


Peck and Hazlewood (14): 


X = 1/u 
l/u + 1/2 
p = 1/۸ + 1/u 


` 1/۸ + 1/u + 


F represents the overall efficiency of the repair system. Then based on 


these values for F and X, the following equations are developed: 


ES 
ዘ 


N [1-F (1- X)] 


= 
II 


(1 - F)/F X 
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To determine MA in terms of E an additional value, D, delay probability; 
is defined by Peck and Hazlewood, and must also be investigated for use 

in solving the queuing problem. The delay probability is the probability 
of an arriving unit not being able to proceed directly into service upon 
arrival at a facility. To obtain this quantity it must be noticed that 

it is no longer true that an arrival is equally likely to occur in any 
interval of time At. The probability of an arrival during At is dependent 


of the number of units in the system. Notice that: 
P(unit arrives for service during At) 


N 
= ) P(unit arrives for service during At/n in the system) En 
n 


(N - n) P. À At 


ll 
Y maz 


Now, let delay probability be defined as follows: 


D = P(unit arrives for service during At and has to wait) 


+ P(unit arrives for service during At) 


N 
) (N - n) P_ A At 
thus, D = Td 
3 GG ee ce 


We then can obtain the probability mass function for the number of units 
in the system ahead of the arriving unit or customer by noticing that 


P(arrival/n in the system)P_ 


P (n in the system/arriving unit) = P (arrival) 


ፇ 
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A (N FE n) P 
or P = E 
P 
(N-n) P 
n=0 x 


In designing this system the service is assumed to be on a first 
come first served basis. Therefore the density function for the waiting 
time can be expressed as 

》 ^ 
f(t) = f(t /n) P 
q 5 6 q n 
where: £ (t, /n) = "089900 EOE خو‎ 0,1,20. 6 and £ (t, /n) is the prob- 
ability distribution of the sum n-c+1 independent and identically dis- 
tributed random variables with the parameters cu where n=c,ctl,...,N-l. 


Thus, £ (t /n) is a gamma distribution with parameter n-ctl and cy if n > 


c or 
c-r 
P E 
=0 n q 
f(t )= 
(የዕ) 
NEL 1] SCHC HE 
| y S C fo + 8001110 
1.4 [1መሙር- ۹ 1 
N ^ 
Since E(t.) = እ ሽ (c.n) P. 


=C 


we see that for n> c 





NS 


n-c 


n-c 
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n-c+l (N-n) 7 


CU N 
ነ (N-n) 


n=0 


n-c+1 
cu 





n! 
(N-n) 210-1 (cp) 


(1/u + 1/۸) 


(1/ህ ተ 1/1 + WË 


N! 
(N-n) (N-n)! 


۲ 
n 


C 


el 


n 


n 


p 


N 
n N! 
十 7 (N-n) KR 


1/u 


(l/u + 1/1 


ር 
ር 
ር! 


| 


0 


n 





APPENDIX B 


B.1 Fortran Program Listing of System Design Model 


B.2 Input Data for System Design Problem 


B.3 Output Results for System Design Problem 
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NNAGE 
NAGE 

D 

XN 

XNC 
XND 
NFCOST 
ALEOC 
ALEOD 
ENDXN 
AGEINC 
PERCNT 
SUPPC 
TRNPC 
Ẹ 

SV 

NP 
ZENI 
PFAC 
FFAC 
LFS 
DYCOST 
DDCOST 
CC 

DC 

ZN 
TUAES 
FAECN 
FAECDN 
WCNT 


TTAENC = 
TTAEND = 


TTAEST 
ھ1‎ 741 
STRN 

RSUPN 
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LIST OF KEY VARIABLES 


= maximum life expectancy of a unit 


current age of a unit for calculation 


= demand 
= total number of units 


number of units in IRL 


= number of units in DRL 


number of cost areas for facility calculations 
first cost of establishing IRL 
first cost if establishing DRL 


= population size for termination of calculations 


vn 


increment of NAGE 
allocation of failed units to IRL and DRL 
percent allocated for recurring support 


= percent allocated for training and init. support 


Initial unit cost 


= salvage value of unit 
= life expectancy of unit 


interest (discount) rate 


= initial cost of IRL facility (PFACD for DRL) 


salvage value of IRL facility (FFACD for DRL) 


= life expectancy for IRL (LFSD for DRL) 


admin, labor, etc. for IRL 
admin., labor, etc. for DRL 
number of repair channels in IRL 
number of repair channels in DRL 


= daily cost of not meeting a demand 


TTAECN = 


total unit annual equivalent cost 
facility annual equivalent cost, IRL 
facility annual equivalent cost, DRL 
yearly cost of not meeting a demand 
total cost of units allocated to IRL 
total cost of units allocated to DRL 
cost of units for total system 

cost of facilities for total system 
training and initial spares cost 
recurring support cost 

grand total annual equivalent cost 
(program search based on this term) 
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SELECTED SENSITIVITY AND CONSTRAINT STUDIES FOR A FINITE 
POPULATION QUEUING SYSTEM DEPLOYED TO MEET A DEMAND 
by 


Carl K. Hinger 


(ABSTRACT) 


The purpose of this thesis is to study the sensitivity of certain 
parameters in a conceptual system containing a finite population of 
units deployed to meet a demand. Selected constraints on the design 
variables are also treated. The measure of effectiveness in this system 
design is the total annual equivalent cost. This cost is the summation 
of the total population unit costs, depot and intermediate repair fac- 
ility costs, training and initial support costs, recurring support costs, 
and cost of not meeting a demand on the system. 

Three separate areas are addressed in this study. The establishment 
of a standard system design, the study of the variation in the cost as 
compared to the standard design, and the effects of constraints on the 
system cost compared to the cost of the standard design. 

A special computer algorithm was developed to perform a multiple 
search through the design variables. This search yields the minimum 
total system cost and a resultant system configuration to be specified 
by the designer which is optimum. The algorithm was then used to analyze 
the effects on the minimum total system cost of variation in system 


parameters and constraints on design variables. 
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